Phase separation in fermionic systems with particle-hole asymmetry 
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We determine the ground-state phase-diagram of a Hubbard Hamiltonian with correlated hopping, 
which is asymmetric under particle-hole transform. By lowering the repulsive Coulomb interaction 
U at appropriate filling and interaction parameters, the ground state separates into a hole and 
an electron conducting phases: two different wave vectors characterize the system and charge- 
charge correlations become incommensurate. By further decreasing U another transition occurs at 
which the hole conducting region becomes insulating, and conventional phase separation takes place. 
Finally, for negative U the whole system eventually becomes a paired insulator. It is speculated 
that such behavior could be at the origin of the incommensurate superconducting phase recently 
discovered in the ID Hirsch model. The exact phase boundaries are calculated in one dimension. 

PACS numbers: 05.30.Fk, 71.10.Fd 
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One very debated issue in the context of correlated 
fermionic materials is the occurrence of phase separation 
(PS). For instance, in ultracold fermionic atoms trapped 
on optical lattices PS appears when the superfluid phase 
coexists with the normal phase P, Q ; these systems may 
be investigated by means of the negative U Hubbard 
model Also, in high-T c materials PS leads to the for- 
mation of one-dimensional stripes It is believed that 
the phenomenon could be a consequence of the interplay 
between antiferromagnetic order and holes propagation; 
as such, its appearance is investigated within the con- 
text of the t — J model [H, @, 0, Hf , which describes the 
Hubbard model in the limit of strong repulsive Coulomb 
interaction. In the present paper we suggest a further sce- 
nario for the occurrence of PS in Hubbard-like systems, 
resulting from competition between conducting phases 



with different charge: in the hole rich region the carri- 
ers are the electrons, while in the pair rich region they 
eventually become the holes. In absence of particle-hole 
symmetry, the proposed mechanism can generate incom- 
mensurability. The mechanism is proven to survive in 
dimension greater than 1 for a simple case. 

We study the bond-charge extended Hubbard Hamilto- 
nian, which has been introduced to describe compounds 
containing extended orbitals Q . It comes from assuming 
that the charge in the bond affects the effective potential 
acting on valence electrons; the extension of the Wan- 
nier orbitals and the hopping between them should vary 
with the charge. A generalization of that Hamiltonian is 
the U — X — X' Hamiltonian introduced in [l(| . In the 
grand-canonical ensemble it reads: 
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where the lower case symbols denote that the coefficients 
of the interactions have been normalized in units of the 
hopping amplitude, and the 3-body interaction term x' 
appears as an effective interaction from the three band- 
model. Moreover N is the number of electrons on the 
-D-dimensional L sites lattice A^. The model belongs to 
the class of Hubbard Hamiltonians with correlated hop- 
ping, which contains some integrable cases that have been 
widely studied in recent years (see ll| and references 
therein) . 

Hbc is in general is not invariant under particle-hole 
transform. For this reason the model has been first pro- 
posed (for x' = 0) in two dimensions, motivated by a 
theory of hole superconductivity [l2| . It was investigated 
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by means of bosonization approach in [13[ at half-filling 
for arbitrary (and weak) x and x' . Also, a particle-hole 
invariant subcase of Hbci which takes place at x' = 2x, 
has been previously discusse dll4ll . 

Recently, it has been proven |l5|.[l6j that in ID at half- 
filling Hbc displays a rich phase diagram already at 
x' = 0, quite different from that of the Hubbard model, 
and from that derived in the weak coupling limit [l3j |. 
In particular, at sufficiently low u > , and for appro- 
priate range of x values an unexpected incommensurate 
superconducting (ICSS) phase appears, which drives the 
usual transition to an insulating state to higher values 
of u. Here we investigate the origin of such phase by 
constructing the exact ground-state phase-diagram of a 
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particular sub-case of (1) still asymmetric under particle- 
hole transform. 

Hamiltonian |T]) can be written in terms of the Hub- 
bard projectors, = \u)i{(3\i\ here \a)i are the states 
allowed at a given site i, a = 0, |, I, 2 (|2) = | f J.)). In- 
serting the choice a: = 1 in Hbc, the resulting Hamilto- 
nian if = Hbc( x = 1) turns out to preserve the number 
of doubly occupied states N d =< N d >=< X),--X"? 2 >: 

[H, Nd] = for arbitrary x', D and N. In the following, 
we shall adopt such choice. H can then be written as 



H = H { 
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where 



Hoi = - E (*r°*f + h.c.) + n e *f (1) 
H 12 = -t x E (A^Af + fc.c.) + ([/ - M ) E *f (?) 

with t x = 1 — x'. We shall limit our analysis to the 
range < t x < 1. By implementing the transformation 
Cj.cr — ► cj^, it can be realized that the range is represen- 
tative of the behavior of the model at any t x value. In 
the two limits t x = (x' = 1) and t x = 1 [x' = 0) the 
Hamiltonian reduces to known cases, namely the infinite 
U Hubbard model, and the bond-charge Hubbard model 
at x = 1 (l7j . both integrablc in one dimension. 
On general grounds, we expect that at large enough posi- 
tive u the ground state would contain the minimum num- 
ber of doubly occupied sites. Hence it will coincide with 
that of the infinite U Hubbard model (no doubly occu- 
pied sites) for N < L, and with its particle-hole coun- 
terpart for N > L. In Fig. [T] such state is denoted as 
Uoo. Also, for large enough negative u, the number of 
doubly occupied sites will be maximized, so that for even 
number of electrons the hopping term will be ineffective 
and the ground state would be a highly degenerate insu- 
lator consisting of N/2 pairs of electrons. In Fig. [T]this 
is denoted as PI (paired insulator). In order to investi- 
gate what happens between these two limits, let us think 
of the Hilbert space as factorized into three orthogonal 
subspaces, 



H = Hqi © H12 © H± 



(3) 



states with no doubly occupied sites belong to Hoi , states 
with no empty sites live in H12, and the remaining states 
stay in Hj_. H10 does not act on doubly occupied sites 
present in H±, and annihilates states belonging to H12', 
hence, at given N d and N, it reaches its minimum in Hoi- 
Analogously, H12 does not act on empty sites present in 
Hj_, and annihilates states in Hoi, hence reaching its 
minimum in Hi2- So that, at any given N and Nd, the 
absolute minima of both H10 and H12 (and consequently 
of H) are reached in the space orthogonal to H±. We 



can then rewrite the ground state energy E gs in the form 

E gs = min WeH (ip\H\ip) = min W£noiQni2 (^\H\i)} , 

(4) 

where the minimum has to be taken with respect to Nd, 
the constraint L — N = N e — Nd is implemented through 
the chemical potential, and N e = ^ . is the (con- 
served) number of empty sites. More explicitly, one could 
look for a ground state of the form: 



mN e ,N d )) gs - a\i> m {N' e )) + y/l-a*\iki{N' d )) , (5) 

where a is a further variational parameter, and N e = 
a 2 N' e , N d - (1 - a 2 )N' d . Here |Voi(A^)> e Hoi, and 
|^/'i2(A r ^)} £ H12 are the ground states of Hoi at given 
N' e and Hi 2 at given N' d respectively. The corresponding 
energies are recovered from those of the infinite U Hub- 
bard model, £^?. In fact, the spectra of Hoi and H12 
are the same: up to a multiplicative factor t x , and to 
the conserved quantities N ei Nd, they coincide with the 
spectrum of the infinite U Hubbard model. Explicitly, 



E gs = a 2 E™(K) + (1 - a 2 ) [t x E™(N d ) + uN' d ] 
— n (Nd — N e + L — N) . 
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By minimization of E gs with respect to /1 one obtains 
that the actual value of a 2 is given by 



a 2 = 



L-N + N' d 



NL 



(7) 



The case a 2 = 1 (a 2 = 0) correspond to the Uoo phase 
described above for N < L (N > L), with N' e = L - N 
(N' d = N — L). Apart from these limits, the system is 
always separated into two phases, characterized by N' e , 
N' d , and fi values which minimize E gs as given by ©, at 
given t x and u. Such values do not depend on the actual 
filling N: PS implies that the chemical potential is con- 
stant with the filling. 

Since t x < 1, independently of the actual value of £™? 
(and D) one could expect three different behavior to 
emerge from the minimization equations within the PS 
phase, depending on the filling and the interaction pa- 
rameters. 

(i) N' e < N' d < L: both the coexisting phases are con- 
ducting. They arc characterized by different Fermi 



momenta kp^ and hp 1 ' , which can be inferred from 
those of the corresponding infinite U model, and are 
in general not commensurate between each other. 
For this reason we characterized the phase as IPS, 
incommensurate phase separation. Later on we 
shall discuss how such incommensurability reflects 
onto some physical feature. 

(ii) N' e < N' d — L; only the hole- rich region is conduct- 
ing, whereas the other phase is a paired insulator. 
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Figure 1: (Color online) Phase diagram in the u-t x plane in 
D = 1: the transition to the Uoo region depends on the fill- 
ing n, whereas the other transitions are independent of it. 
CPS and IPS denote the conventional and the incommensu- 
rate phase separation regions respectively; PI stays for paired 
insulator. Full dots indicate the tricritical points. 



The system is characterized by the Fermi momen- 
tum kp^ of the conducting electrons, i.e. that of the 
infinite U Hubbard model at a filling L — N' e < N . 
We denote such phase as CPS (conventional phase 
separation) . 

(iii) N' e = N' d = L; both phases arc insulating. In 
this case we are in one of the infinitely many (in 
the thermodynamic limit) degenerate ground states 
characterizing the PI phase. 

A further interesting limiting case, not discussed above 
though contained in the minimization equations, occurs 
when N' e = N' d ^ L. In such case both phases are con- 
ducting as in (i), but they are characterized by the same 
Fermi momentum, kp = kp . This is the case for t x = 1. 
The latter is the already mentioned bond-charge Hub- 
bard model at x = 1 [121], and manifests particle-hole 
invariance. This example illustrates how in our scheme 
the IPS phase, intended as the appearance of two not 
commensurate modulations, is in fact associated with the 
presence of particle-hole asymmetry. 

In order to make the above predictions more quanti- 
tative, we now turn to the discussion of the ID case. In 
that case the ground-state energy per-sitc of the infinite u 
Hubbard model is known exactly. In the thermodynamic 
limit it reads: 



■ sin (jm) 



(8) 



with n average per-site filling, n = 

When inserted into E gs /L, as given by ©, the above in- 
formation allows for an explicit derivation of the ground- 
state phase diagram, shown in Fig. [T] This is achieved 
by solving the minimization equations with respect to 



JV' 



and n'j 



L 



which in turn fix k 



(e) 



and kp = 7r(l — n' d ). 

Interestingly, the same phase diagram and ground-state 
energy can be also obtained by generalizing the approach 
developed in [l8| for studying PS in the extended Hub- 
bard model for V — > oo, with V neighboring sites diag- 
onal Coulomb interaction. At a given filling, the ground 
state can be thought of as an island of Nd + Z singly and 
doubly occupied sites inserted into a chain of L — Nd — Z 
singly and empty sites. Such choice allows to maximize 
the number of available low energy momenta. In the ther- 
modynamic limit, minimizing the energy of that state 
with respect to Nd and Z gives precisely the same phase 
diagram obtained here in Fig. 1. It turns out that the 
following relations hold: n' d = N £ z , and n' e = j^w+z • 
The solutions n' e = 1 — n (a 2 = 1) below half- filling, 
and n' d = 1 (cr 



0) above half-filling, are obtained for 
u > u c (n,t x ) (see below). These identify the Uoo phase. 
In particular, at n = 1 a charge gap A c opens and the 
phase becomes insulating: for u > u c (\,t x ) — 2(t x + 1), 
A c = u — u c . 

For a 2 ^ 0,1 the solution to the minimization equation 
reads: 



, 1 

n„ = — arccos 

7T 



(f) 



- arccos ( 
1 



IPS 
CPS 



(9) 

Minimizing E gs is thus reduced to solving the following 
transcendental equation in ji, 



1 



/' 



— (—2t x sin irn'd + un' d + 2 sin im' e ) 



(10) 



with n' e and n' d as specified in ©. Let us denote the 
solution by p,(u,t x ). The actual value of /2/ in the IPS 
phases is limited by the constraint 



Hi G ]max(— 2, u — 2t x ) 7 min(2, u + 2t x )[ 



(11) 



Whereas in the CPS phase fj, becomes independent of t x , 
and p,c(u) is the solution of (flO|) for n' d = 1. In fact, all 
the transition curves can be characterized by the value 
of p, along them: 



1. Uoo -> IPS, CPS, 



fj,(u,t x ) 



— 2 COs(7Tn) 77 < 1 

u + 2t x cos(7m) n > 1 



(12) 



7T(l-<) 



2. IPS — > CPS, [ii = uc] the transition occurs at 
tx = \fic(u) —u]/2, implying n' d = 1. 

3. finally at CPS -> PI JL C = -2. 

When passing from fj, to n, one realizes that only the 
two transitions 1 do depend on the filling, whereas the 
others are independent of it. At transition 1 the ex- 
plicit value of u c (n,t x ), is obtained by solving eq. (fT2")) 
for u. Transition 2 takes place at t x = i x ; transition 
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3 occurs at u = —4. As an example, we report the 
critical curves in fig. [1] in the (u,t x ) plane at two dif- 
ferent filling values. For u < 2, depending on n, two 
tricritical point (Tl) are recognized, characterized by 
the merging of Uoo, IPS and CPS phases: one (not 
shown) occurs at n = 2, whereas the other one occurs 
at the same u value and at a filling npi < 1 . Ana- 
lytical calculations show that Utx = — 2(cos[7rnri] + t x ), 
n-ri ~ \ arccos[(3t2; — 1)/(1 + t x )]. Increasing t x pro- 
gressively drives the system to a particle-hole symmetric 
diagram. For t x = 1 npi — 0, which is the particlc-holc 
counterpart of the tricritical point in n = 2. The criti- 
cal curves in the (u,ri) plane (not shown) highlight that 
particle- hole asymmetry (tx ^ 1) favors the IPS phase at 
filling greater than half, in which case it survives at ap- 
propriate positive u values up to n = 2. Fig. [T]also shows 
that the infinite U Hubbard model (t x = 0) exhibits CPS 
for u < 2 if doubly occupied sites are allowed. 

In order to clarify the meaning of the different mod- 
ulations in the IPS phase, we also give in Fig. [5] the 
charge-charge correlations C(r) = ((n, — n)(nj+ r — n)) 
at half filling (upper part), and their Fourier transform 
N(q) J2 r eiqr C(r) (lower part). Since C(r) can be eval- 
uated exactly (see for instance we can also provide 
an analytic expression for N(q), which reads: 



N(q) = C(0) + (5(q)-l)d 



ps-a 2 7 fe (") {q)-( l - a2 )i k w (q) 

F (13) 

Here S(q) is the Dirac delta, and dps = (1 — n + n' d )(n — 
1 + n' e ) is the constant responsible for the divergent con- 
tribution in q = characteristic of phase separation: at 
given (u, t x ) it becomes zero for n < ni = 1 — n' e , or 
n > rih = 1 + n'j. Moreover 



lk{q) 



4tt 



;K[L 2 (e 



*( 9 +2fc)) + L2 ( el ( g -2fc))_ 2j L 2 (e i;<? )] , 

(14) 

with k = kp,kp, and L 2 (z) denoting the dilogarithm 
of z (R[] being the real part of []), describes the q- 
dependence of density correlations. Equations (fT3")) - (fl"l|) 
enlighten how in the IPS phase two wave- vectors do char- 
acterize N(q), namely q^ = 2kp^ and q^ = 2k^p\ as 
shown in the lower part of Fig. O Whereas in the CPS 
phase, since k^ = and 70(9) = 0, only one feature 
characterizes N(q), namely the wave- vector of the single 
fermion problem corresponding to a filling n = 1 — n' e . 
It can be recognized that the above type of ground-state 
phase diagram could explain different features displayed 
by the non-integrable ID Hirsch model discussed in [lq . 
At variance with the present model, in that case the 
model still had arbitrary t x (t x = 2x— 1) but also x' = 0, 
implying the appearance of a further term in the Hamil- 
tonian. The numerical solution did show incommensu- 
rate charge-charge correlations (and superconductivity) 
at appropriate values of u and t x , always missing particle- 
hole invariance. Also, by suitably varying u and t x at 
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Figure 2: (Color online) (a) Equal time density correlations 
C(r) vs distance r at u = 1 and t x — 0.6; the dashed line 
represents the slope 1/r 2 . (b) N(q) — dps at t x = .3, u = 0, 
and three different fillings: from top to bottom n = .5 (blue), 
n = 1 (pink), n — 1.7 (violet). 



half-filling, the system could first pass from an insulat- 
ing phase, characterized by flat density correlations, to 
a conducting phase in which that correlations acquire 
a modulation; then to the ICSS phase in which further 
modulations occur. We suggest that the two latter phases 
could be identified with our CPS and IPS phases respec- 
tively. As an example, in fig. [2] we report C(r) at u and 
t x values identical with those used for the Hirsch model 
in fig. 3 of The resemblance of the two figures is 

convincing. 

In this paper, we have shown how the appearance of 
incommensurate characteristic wave-vectors in strongly 
interacting fermionic systems could be related to the oc- 
currence of phase separation, in case the system is not 
symmetric under particlc-holc transform. Recently the 
occurrence of incommensurate modulations in charge cor- 
relations was observed also in a Hubbard-like model in- 
cluding next-nearest neighbors hopping (20|. Since that 
model is particle-hole asymmetric, we suggest that IPS 
could occur in that case. 

The mechanism could also be a natural candidate for 
explaining PS in unbalanced mixtures of cold fermionic 
atoms[H [ij]. In this case the PI phase should be in- 
terpreted as the (unpolarized) superfluid phase, and the 
CPS phase should be identified with coexistence of the 
superfluid with the (polarized) Fermi gas. 
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